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EXECUTIVE  SUMMARY 


OBJECTIVE 

This  report  proposes  a  discrete-time  model  of  an  Asynchronous  Transfer  Mode  (ATM)  statistical 
multiplexer  subject  to  heterogeneous  groups  of  on-off  traffic  sources  with  correlated  burst  periods. 
NRaD  engineers  have  formulated  algorithms  to  compute  both  the  exact  transient  and  steady-state 
cell-loss  probabilities  of  the  multiplexer. 

METHOD 

The  complexity  of  the  algorithms  is  exponential  to  the  number  of  traffic  types  multiplexed;  there¬ 
fore,  it  may  be  impractical  to  use  them  to  address  realistic  ATM  design  problems.  For  this  reason, 
NRaD  engineers  formulated  and  analyzed  an  approximation  model  of  the  multiplexer  subject  on-off 
traffic  sources  with  random  burst  periods.  They  developed  an  algorithm  to  calculate  the  exact  cell- 
loss  probability  of  the  approximation  multiplexer  model.  Unfortunately,  the  time  complexity  of  the 
algorithm  is  also  exponential  to  the  number  of  traffic  types  multiplexed.  To  circumvent  this  difficulty, 
the  engineers  formulated  an  efficient  approximation  procedure  to  estimate  the  cell-loss  probability  of 
the  approximation  multiplexer  model. 

CONCLUSION 

NRaD  engineers  developed  algorithms  to  compute  both  multiplexer  models’  cell-loss  probabilities. 
The  algorithms’  results  have  been  compared  to  simulation  outputs  and  they  have  showed  perfect 
agreement  (within  noise  level).  The  algorithms  will  be  invaluable  for  resolving  basic  ATM  design 
issues  such  as  buffer  dimensioning,  maximal  loading,  and  admission  control. 
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1.  BACKGROUND 


Asynchronous  Transfer  Mode  (ATM)  networks  will  dominate  future  global  telecommunication 
industries.  These  networks  will  make  many  telecommunication  services  (e.g.,  video  conferencing  and 
high-speed  data)  economical  and,  in  turn,  will  bring  fundamental  changes  to  the  way  that  people  com¬ 
municate  and  travel.  Despite  the  expected  prominence  of  ATM,  many  basic  ATM  design  issues,  such 
as  buffer  dimensioning  and  routing,  still  have  not  been  addressed.  The  resolution  of  these  issues  will 
require  a  thorough  understanding  of  the  dynamic  behavior  of  the  ATM  multiplexer.  The  analysis  of 
the  cell  loss  behavior  of  an  ATM  multiplexer  is  the  focus  of  this  report. 

2.  MATHEMATICAL  MODEL  OF  THE  ATM  MULTIPLEXER 


One  classical  approach  to  understanding  a  system’s  behavior  is  to  design  a  mathematical  model  that 
captures  the  salient  characteristics  of  the  system  and  then  analyze  the  model.  This  approach  will  be 
used  in  this  report. 

Statistical  multiplexing  is  a  means  of  sharing  a  communication  channel  by  many  traffic  sources 
simultaneously.  A  statistical  multiplexer  consists  of  a  channel,  a  buffer,  and  a  set  of  traffic  sources 
sharing  the  channel.  Traffic  is  generated  by  the  sources  and  removed  by  the  channel.  Whenever  the 
instantaneous  aggregate  traffic  rate  of  the  traffic  sources  is  below  the  channel  capacity,  the  amount  of 
the  traffic  in  the  buffer  does  not  accumulate;  however,  when  the  total-traffic  rate  exceeds  the  capacity 
of  the  channel,  the  excess  traffic  will  be  stored  in  the  buffer;  and  when  the  buffer  is  full,  excess  traffic 
will  be  lost. 

The  most  difficult  task  in  developing  an  accurate  ATM  statistical  multiplexer  model  is  the  charac¬ 
terization  of  a  traffic  source.  As  one  may  expect,  the  traffic  generation  behavior  of  a  traffic  source  is 
stochastic  and  differs  depending  on  the  traffic  type.  The  problem  is  to  develop  a  general  traffic  model 
that  may  be  used  to  describe  an  arbitrary  traffic  source.  The  graph  (figure  1)  depicting  the  traffic  gen¬ 
eration  rate  of  a  traffic  source  as  a  function  of  time  yields  insight  for  the  formulation  of  a  satisfactory 
traffic-source  model.  Common  to  all  bursty  traffic  sources,  there  are  burst  periods  during  which  a 
traffic  source  would  generate  traffic  close  to  some  peak  rate,  P,  and  there  are  idle  periods  during 
which  the  source  would  generate  close  to  nothing.  While  it  is  impossible  to  completely  describe  all 
the  particular  characteristics  of  a  traffic  source,  there  are  three  fundamental  characteristics  that  are 
considered  to  be  the  most  important  in  determining  the  cell  loss  behavior  at  a  statistical  multiplexer: 
the  peak  traffic  rate,  the  mean  traffic  rate,  and  the  expected  burst  period.  It  is  imperative  that  a  traffic 
source  model  would  account  for  these  three  characteristics. 

An  arbitrary  traffic  source  that  has  a  peak  traffic  rate,  P,  a  mean  traffic  rate,  m,  and  an  expected 
burst  period,  E(BP),  may  be  approximated  by  an  on-off  traffic  source  model  that  is  described  as  fol- 
lows: 

An  on-off  traffic  source  of  the  i-th  type  alternates  between  on  and  off  periods.  The  duration  of  an 
on  period,  measured  in  terms  of  integral  multiple  of  a  time  slot  of  duration  r,  is  geometrically  distrib¬ 
uted  with  a  parameter  of  a,;  the  duration  of  an  off  period  is  geometrically  distributed  with  a  parameter 
°f  A,  where  0  <  o., ,  /?.  <  1 .  During  its  on  periods,  the  source  generates  Px  cells  per  time  slot  and  dur- 
ing  its  off  periods,  generates  nothing.  A  schematic  of  a  traffic  source  of  the  i-th  type  is  shown  in  fig¬ 
ure  2. 
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Figure  1 .  Traffic  generation  rate  of  a  traffic  source. 
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Figure  2.  Traffic  generation  rate  of  an  i-th  type  on-off  source. 

Notice  that,  given  an  arbitrary  traffic  source  with  peak  rate,  P,  mean  rate,  m,  and  expected  burst 
period  E(BP),  one  can  construct  an  on-off  traffic  source  approximation  that  has  the  same  three  funda¬ 
mental  characteristics.  The  parameters  P(  ,  a,-,  and  of  the  approximating  on-off  source  are  uniquely 
determined  according  to  the  following  formulae: 


*  a'  =  1-£^) 

*  Pi  =  1 - l~p - 

Because  of  its  extreme  nature  (generating  at  peak  rate  during  its  on  period  and  nothing  during  its 
off  periods),  an  on-off  traffic  source  with  fundamental  characteristics  P,  m,  and  E(BP)  is  likely  to  be 
more  bursty  than  most  traffic  sources  with  the  same  characteristics;  therefore,  it  would  effect  a  higher 
cell  loss  at  a  multiplexer.  This  is  particularly  desirable  for  modeling  purposes. 
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3.  TRANSIENT  CELL-LOSS  PROBABILITY 

A  traffic  type  typically  has  two  requirements:  delay  and  loss.  The  delay  requirement  of  a  traffic 
type  is  expected  to  be  met  by  an  ATM  network.  The  loss  requirement  is  a  bigger  concern  for  the  net¬ 
work.  For  this  reason,  the  loss  characteristics  are  the  most  important  performance  measure  of  an  ATM 
multiplexer.  The  focus  of  this  work  is  to  formulate  an  efficient  procedure  that  computes  the  transient, 
as  well  as  steady-state,  cell-loss  probability  of  the  multiplexer.  A  transient  solution  is  desirable  since 
it  is  expected  that  the  behavior  of  an  ATM  network  is  highly  dynamic  and  a  steady-state  solution 
alone  may  not  provide  adequate  information  for  the  control  and  design  of  an  ATM  network. 


4.  LITERATURE  REVIEW 

Statistical  multiplexing,  an  efficient  and  flexible  means  of  sharing  a  communication  channel  by 
many  traffic  sources  simultaneously,  is  so  critical  to  an  ATM  network  that  the  term  asynchronous 
transfer  mode  is  just  an  alias  for  statistical  multiplexing.  Because  a  thorough  understanding  of  the 
effect  of  statistically  multiplexing  diverse  traffic  sources  onto  a  channel  is  required  to  address  many 
basic  ATM  design  problems,  numerous  attempts  were  made  to  analyze  the  multiplexer’s  behavior. 
The  continuous-time  analog  of  the  on-off  model  developed  in  this  report  was  first  proposed  in  Anick, 
Mitra,  and  Sondhi  (1982).  The  main  motivation  of  that  pioneering  work  was  to  evaluate  the  steady- 
state  cell  loss  behavior  of  a  statistical  multiplexer  subject  to  only  one  type  of  traffic  stream.  A  closed- 
form  expression  for  the  steady-state  distribution  of  the  buffer  content  was  derived  assuming  that  the 
buffer  size  of  the  multiplexer  is  infinite.  Although  the  analysis  is  elegant,  it  cannot  be  used  in  the  con¬ 
trol  and  design  of  an  ATM  network  since  only  one  traffic  type  was  allowed  in  that  model,  whereas  in 
a  realistic  ATM  network,  heterogeneous  traffic  types  are  expected. 

In  Kosten  (1984),  the  on-off  traffic  model  is  extended  to  include  heterogeneous  types  of  on-off 
traffic  sources.  An  approximation  expression  to  estimate  the  steady-state  probability  of  the  buffer 
content  exceeding  a  very  large  value  is  derived  under  the  assumption  that  the  buffer  size  is  infinite. 
The  validity  of  his  approximation  had  never  been  extensively  verified  by  simulation  studies. 

Kobayashi  and  Ren  (1992)  investigated  a  multiplexer  subject  to  the  same  assumptions  in  Kosten 
(1984),  but  they  focused  on  obtaining  the  transient  probability  of  the  buffer  content  exceeding  a  cer¬ 
tain  value.  Unfortunately,  they  gave  only  a  formulation,  rather  than  an  explicit  solution  of  the  prob¬ 
lem.  Shum  (1994)  derived  an  approximation  expression  for  transient  probability  of  the  buffer  content 
exceeding  a  large  value;  however,  the  accuracy  of  his  expression  is  not  uniformly  satisfactory.  To 
obtain  an  exact  solution  of  the  distribution  of  buffer  content  exceeding  a  certain  value  for  the  multi¬ 
plexer  model  with  an  infinite  size  buffer  requires  one  to  solve  a  set  of  partial  differential  equations 
with  very  specialized  boundary  conditions;  unfortunately,  not  even  a  satisfactory  numerical  procedure 
exists  to  solve  these  equations.  A  model  with  &  finite  buffer  is  even  less  tractable. 

Currently,  a  satisfactory  procedure  that  estimates  the  cell  loss  rate  effected  by  heterogeneous  on-off 
traffic  sources  at  a  finite-buffer  statistical  multiplexer  does  not  exist,  and  developing  such  a  procedure 
is  the  main  motivation  for  this  report. 
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5.  ANALYSIS 


We  define  the  following  variables: 

•  r  the  duration  of  a  time  slot 

•  M  the  number  of  different  traffic  types 

•  Nj  the  number  of  traffic  sources  of  the  i-th  type  i= 1 , ...,  M 

•  R[  the  number  of  cells  generated  by  an  i-th  type  traffic  source  in  a  time  slot  during  the  source’s 
on  periods 

•  B  the  maximum  number  of  cells  that  may  be  buffered  at  the  multiplexer 

•  C  the  number  of  cells  that  may  be  transmitted  by  the  output  link  of  the  multiplexer  per  time 
slot;  We  assume  1  <  C 


We  make  the  following  simplifying  assumptions: 

a.  An  on  period  of  a  traffic  source  the  i-th  type  lasts  Xt  time  slots,  where  Pr{Xt  =  j )  =  a^O  -  a,  ) 
and  1  <  j  and  0  <  a,  <  1 

b.  An  off  period  of  a  traffic  source  the  i-th  type  lasts  Yt  time  slots,  where  Pr(Yt  =  j )  =  fa^il  -fa) 
and  1  <  j  and  0  <  fa  <  1 

c.  The  number  of  cells  at  the  buffer  at  a  time  k+1  is  computed  as  max(min(z  +R  ■  J  -  C,  B), 0), 
where  i  is  the  number  of  cells  in  the  buffer  at  time  k;  J  is  the  state  of  the  traffic  sources  at  the 
start  of  time  k 


m  .j  _ 

d.  (Stationary  Condition)  Y '  — - — ~a'  —  <  c 

*  J  A  i  1 

l-a,+  l-A 

M 

e.  >  C 

i  —  1 

Assumptions  a.  and  b.  state  that  the  on  and  off  periods  of  an  i-th  type  traffic  source  are  geometri¬ 
cally  distributed  with  parameters  at  and  fa,  respectively.  Assumption  c.  states  that  the  number  of  cells 
at  the  beginning  of  the  (k+l)-th  time  slot  is  equal  to  i,  the  number  of  cells  at  the  buffer  at  the  begin¬ 
ning  of  the  k-th  time  slot,  plus  the  number  of  cells  generated  during  the  k-th  slot,  R  ■  J,  subtract  the 
number  of  cells  removed  by  the  channel  during  the  time  slot,  C;  however,  the  number  of  cells  in  the 
buffer  will  never  exceed  the  buffer  size  B  and  it  will  never  drop  below  0.  Assumption  d.  states  that 
the  overall  mean  aggregate  traffic  rate  of  all  sources  is  less  than  the  capacity;  this  allows  the  number 
of  cells  at  the  multiplexer  to  have  a  steady-state  distribution.  Assumption  e.  ensures  that  we  do  not 

M 

have  a  trivial  problem,  as  —  C  implies  that  cell  loss  would  never  occur. 


We  denote: 


•  Lk  the  number  of  cells  in  the  buffer  at  the  beginning  of  the  k-th  time  slot 

lim 

•  L  =  k  00  4  the  number  of  cells  in  the  buffer  at  the  beginning  of  a  time  slot 

•  Jik  the  number  of  the  i-th  type  traffic  sources  that  are  on  at  the  beginning  of  k-th  time  slot 
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•  S(M, N)  =  {(/! ,j2, -JM)  1  0  =£  Ji  -  i  =  1, 2,  ...,M\ 

•  S(B,M,  N)  =  j (ij)  I  0  <  i  <  5,)  £  5(M,A0) 

•  Z*  =  (J lb  Jib  ■■■’■! Mk)  ^  S(M,N) 

•  Pr(lk+ 1  =  4+i  1  4  =  4)  ^  -  4  for  4+i  and  4  e  S(M,N)  the  probability  that  7t+1  =  Zjfc+1, 
given  that  Jk  =  lk 

•  (Lk  =  j,Jk  =  /)  the  joint  event  of  Lk  =  j  and  Jk  =  l. 

•  Pr(Lk+  j  =  j,Jk+ 1  =  l  I  Lk  =  =  m)  =  Pk(j,l  I  i,m)  the  probability  that  the  state  at  time 

k+1  is  (/',/),  given  that  the  state  at  time  k  is  (i,m) 

The  transition  probability  from  the  state  (Lk  =  i,Jk  =  m)  at  time  k  to  the  state 
(Lk+ j  =  j,Jk+i  =  /)  at  time  k+1  is 


Pk(j,l  I  i,m)  =  $ 


M 

j  =  max(min(/  +  ^  Rxmx-C,  B),  0) 


*=  1 


Pr{4+i  =  l-'ijc  =  ^}’ 


0) 


where  '£,(A)  is  the  indicator  function  of  the  event  A,  i.e., 


£{A}  = 


1  A  is  true 
0  A  is  not  true 


The  transition  probability  for  the  state  of  the  traffic  sources  from  time  k  to  £+1, 
Pr{Z*+,  =  L  I  Jk  =  m\,  is  derived  as  follows: 


For  now,  we  only  consider  the  change  in  the  number  of  the  on  sources  of  the  h-th  type.  Suppose  at 
time  k,  the  number  of  on  sources  of  the  h-th  type  is  i.  Let  ah  k  denote  the  state  of  a  traffic  source  of 
the  h-type  at  time  k;  that  is,  oh  k  is  either  on  or  off.  We  conclude  that 

Pr(°h.k+i  =  off  I  ah  k  =  on)  =  1  -  ah, 

Pr(°h,k+ i  =  on  I  ohk  =  on)  =  ah, 

Pr(°h,k+i  =  on  I  ahk  =  off)  =  1  -0h, 

Pr(°k,k+ 1  =  off  I  oh  k  =  off)  =  . 

Let  m  be  the  number  of  the  i  on  sources  that  stay  on  at  time  k+1 ;  hence,  m  <  i  and  i  -  m  sources 
turn  off  at  time  k+1.  Let  l  be  the  number  of  the  -  i  off  sources  that  turn  on  at  time  k+1,  where  l  < 
Nh-i.  We  conclude  that  Nh-i-l  of  the  on  sources  at  time  k  stay  off  at  time  k+1 .  Since  the  traffic 
sources  turn  on  and  off  independently,  we  obtain 


where 


Pr(Jh,k+\  ~  j  I  4i,it  ~  0  ~ 

i  E  (mK>*(1  "  7  -  A.)'. 

m  —  \  {U—m—jj '  7  '  / 


x! 

y!(y  -  x)\ 


if  y  <  x  and 


=  0  otherwise. 


(2) 


When  N  is  large,  it  is  preferable  to  use  Stirling’s  formula  to  approximate  N\  by  NNe~Nj2ttN. 
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(3) 


We  conclude  that  since  the  traffic  sources  turn  on  and  off  independent  of  each  other, 

M 

Pr(Ik+l  =  H  Lk  =  m)  =  0^  =  /,  I  Juk  =  m,-) . 

i=i 

If  the  2-tuple  (i,j)  E  S(B,M,N)  were  used  to  denote  as  a  state  of  the  multiplexer,  then  the  set  of 
states  S(B,  M,  N)  forms  a  Markov  chain.  The  transition  probability  from  state  (i,m)  at  time  k  to  state 
(/',/)  at  time  k+1  for  the  Markov  chain  is 

pRLk* i  =  hi *+i  =  L  1  h  =  i,Jk  =  m)  =  £{{/'  =  max(min(7  +  R  ■  m-C, B), OjJPrjj^+j  =  l  I  Jk  =  m}.  (4) 
The  .state  evolution  equation  of  the  multiplexer  is 

pR^t+i  =  j>It+ 1  ~  D  =  I  PRLk+i  =  j’lk+i  ~  D  \  Lk  =  i,Jj.  =  m)Pr{Lk  =  =  m),  (3) 

(i.m)e5(5.A/,A0 

where  L  and  7  denote  the  number  of  cells  in  the  buffer  and  the  state  of  the  traffic  sources  in  steady 
state,  respectively. 

Given  that  the  multiplexer  is  at  state  (i,m)  e  S(B,  M,  N)  at  time  0,  the  probability  that  the  state  of  the 
multiplexer  at  time  k  is  (/,/),  that  is,  the  transient  probability  that  the  multiplexer  is  at  (/,/),  may  be 
recursively  determined  according  to  (4),  if  one  sets  Pr(Ln  =  i,ju  =  m)  =  1  and  Pr(L0  =  x,  J()  =  y)  =0 

M  ~ 

for  (x,y)  *  (i,m).  The  time  complexity  of  obtaining  Pr(Lk  =  j,Jk  =  £)  is  Opfc(fi]^7/,)2). 

/=  1 

The  probability  that  the  state  of  the  multiplexer  is  at  (J,  /)  at  a  random  epoch  sufficiently  far  away 
from  time  0,  that  is,  the  steady-state  probability  that  the  multiplexer  is  at  satisfies  the  following 
fixed-point  equation: 


Pr(L  =j,J  =  0  =  ^  PHL  =  7,7  =  l  \  L  =  i,J  =  m)Pr(L  =  i,J  =  m) ,  (6) 

(i.m)ES(BM.M) 

where  L  and  7  denote  the  number  of  cells  in  the  buffer  and  the  state  of  the  traffic  sources  in  steady 
state,  respectively.  Although  the  probability  Pr(L  =  j,J  =  /)  may  be  computed  recursively  using  equa¬ 
tion  (4)  by  letting  k  oo ,  doing  so  would  take  a  prohibitively  large  number  of  iterations.  The  steady- 
state  probabilities  may  be  obtained  via  the  following  iterative  procedure: 

k  =  0; 

for «i,i)e  S(BMW) 


{ 

^  ~  i  sWJTn]  i  : 

1 ; 

for((i  ,£)  €  S(B  M  JR_)) 
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( 

0>k+%.i)  =  0; 
) ; 


While  ( 

{ 

k  =  k+ 1  ; 

for  ((/,/)  e  S(jBM£)) 


Pk+\ij>-PkU)  -  5  for  each  (itf  e  S(B  M  £)) 


pick  any  e*  such  that  0  <  e*  <  1; 
Pk(i.D  =  (f  ~  £k)pk  1  (i ,2>  +  £* 


X  d  i  -  max(nun(£  +  R-m  -  C £ )0) i J~[/V (7;  =  jt  J-,  =  m^p* 

CfD  1J  \T\  \  f  — 


(es>)  € 


J  i=i 


} ; 


5  =  X  Pk,.L> 

\i&  €  SiJBMJf) 


for  (0',i)  e  S(BMJN)) 


{ 

Pk«.D  = 


PVz). 

5  ’ 


}; 


The  procedure  terminates  whenever  it  has  found  a  consistent  set  of  p^’s  that  satisfy  equation  (6). 


The  condition  ^  RiNi — - — -■  —  - — <  C  ensures  that  the  Markov  chain  is  ergodic.  The  ergodicity  in 

+  I-A 

turn  ensures  that  a  unique  steady-state  distribution  solution  exists.  This  implies  that,  if  the  procedure 
converges,  it  would  converge  to  a  unique  solution.  In  general,  existence  of  a  solution  does  not  imply 
convergence.  Our  computational  experiences  suggest  that  the  procedure  always  converges,  however. 
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M  m 

The  space  complexity  of  the  procedure  is  0(B  J~[  N,);  the  time  complexity  is  0(B2  Nj). 

/=!  /=1 

Once  that  the  state  probability  distribution  of  the  multiplexer  is  determined,  the  cell  loss  rate  at  the 
multiplexer  may  be  found  as  follows: 

We  define  the  following  variables: 

•  NL{L  =  i,J  =  Z)  the  number  of  cells  that  are  lost  at  the  end  of  a  time  slot,  given  that,  at  the  start 
of  the  time  slot,  the  number  of  cells  in  the  buffer  is  i  and  the  state  of  the  traffic  sources  is  / 

•  NA(L  =  i,J  =  D  the  number  of  cells  that  arrive  during  a  time  slot  k,  given  that,  at  the  start  of  the 
time  slot,  the  number  of  cells  in  the  buffer  is  i  and  the  state  of  the  traffic  sources  is  l 

•  LR(L  =  i,J  =  D  the  cell  loss  rate  during  a  time  slot,  given  that  the  number  of  cells  in  the  buffer 
is  i  at  the  start  of  k  and  the  state  of  the  traffic  sources  is  l 

Note  the  following: 


M 

NL(L  =  i,J  =  /)  =  max(^  Rjlr(B-i  +  0,0) 

j-  1 


M 

NA(L  =  i,J  =  /)  =  J^R/j 
/'=! 


LR(L  =  i,J  =  l) 

The  cell-loss  rate  during  a  time  slot  is: 


NL(L  =  i,J  =  /) 
NA(L  =  i,J  =  IJ 


(V) 

(8) 

(9) 


^  LR(L  =  i,J  =  j)  Pr(L  =  i,J  =  f).  (10) 

{ij)eS(BM, N) 

Equations  (7)  to  (10)  outline  the  steps  to  compute  the  steady-state  cell-loss  probability  of  the  multi¬ 
plexer  once  the  steady-state  probability  of  the  multiplexer  at  each  state  (i,m)  E  S(B,M,£D  is  found. 
Equations  analogous  to  (7)  to  (10)  may  be  developed  to  compute  the  transient  cell-loss  probability  of 
the  multiplexer  once  the  transient-state  probability  of  the  multiplexer  at  each  state  (/,m)  e  S(B,M,  AQ  is 
known. 

6.  NUMERICAL  EXAMPLES 


We  conducted  simulation  studies  for  two  arbitrarily  selected  ATM  multiplexing  scenarios  to  evalu¬ 
ate  the  steady-state  cell-loss  probabilities  projected  by  the  algorithm  (see  tables  1  and  2).  The  results 
from  the  simulation  studies  and  the  algorithm  are  as  follows: 

6.1  EXAMPLE  #1 

B  =  2 
C  =  40 
M  =  2 

Nj  =  2  at  =  0.2  /?7  =  0.3  R}  =  20 
N2  =  1  a2  =  0.8  /32  =  0.7  R2  =  10 

We  denote  P(l[)  as  the  steady-state  probability  that  the  buffer  content  of  the  multiplexer  at  the  start 
of  a  time  slot  is  i  and  that  the  state  of  the  traffic  sources  is  j  and  e  as  the  steady-state  cell-loss  proba- 

bihty.  Relative  error  is  defined  as  100  x - - - . 

simulation  result 
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Table  1.  ATM  multiplexing  algorithm  example  #1 . 


Analysis 

Simulation 

Relative  Error  (%) 

P(0,(0,0)) 

0.09659475 

0.09248844 

-4.44 

p(0,(0,1)) 

0.10355416 

0.11048619 

+6.27 

P(0,(1,0)) 

0.19054385 

0.18997625 

-0.30 

P(0.(1,1)) 

0.26508439 

0.26171729 

-1.29 

p(0,(2,0)) 

0.08605149 

0.08351455 

-2.91 

P(0,(2,1)) 

0.12645847 

0.12860892 

+1.67 

P(1,(0,0)) 

0.00000004 

0.00000000 

N/A 

pd,(o,D) 

0.00000002 

0.00000000 

N/A 

P(1,(1,0)) 

0.00000002 

0.00000000 

N/A 

p(1,(1,1)) 

0.00000001 

0.00000000 

N/A 

P(1,(2,0)) 

0.00000000 

0.00000000 

N/A 

P(1,(2,1)) 

0.00000000 

0.00000000 

N/A 

P(2,(0,0)) 

0.01720326 

0.01587302 

-8.38 

P(2.(0,1)) 

0.06709292 

0.07074116 

+5.16 

P(2,(1 ,0)) 

0.00860163 

0.00824897 

-4.28 

P(2,(1,1)) 

0.03354646 

0.03374578 

+0.59 

P(2,(2,0)) 

0.00107521 

0.00024997 

-330.00 

P(2,(2,1)) 

0.00419331 

0.00424947 

+1.33 

£ 

0.042713551 

0.043483995 

+1.77 

6.2  EXAMPLE  #2 

B  =  20 
C  =  27 
M  =  2 

Nj  =2aj  =  0.2  pi  =  0.5  R,  =20 
N2  =  3  a2  -  0.1  p2  =  0.8  R2  =  15 

We  denote  Pi  as  the  steady-state  probability  that  the  buffer  content  of  the  multiplexer  at  the  start  of 

a  time  slot  is  i  and  e  as  the  steady-state  cell-loss  probability.  Relative  error  is  defined  as 

mn  v  simulation  result  -  analysis  result 

1  (JO  X  — — - — — - . 

simulation  result 
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Table  2.  ATM  multiplexing  algorithm  example  #2. 


Analysis 

Simulation 

Relative  Error  (%) 

Po 

0.42601402 

0.425331493 

-0.16 

Pi 

0.05139291 

0.050258995 

-2.26 

p2 

0.00656487 

0.006659867 

+1.43 

p3 

0.01255233 

0.012419752 

-1.07 

Pa 

0.00684444 

0.006959861 

+1 .66 

Ps 

0.00162335 

0.001459971 

-11.19 

Pq 

0.01994044 

0.021119578 

+5.58 

Pi 

0.00242924 

0.002219956 

-9.43 

Ps 

0.11457789 

0.113537729 

-0.92 

PS 

0.01794887 

0.017779644 

-0.95 

Pio 

0.00251679 

0.002519950 

+0.13 

Pii 

0.00524285 

0.005039899 

-4.03 

^12 

0.00245564 

0.002339953 

-4.94 

^13 

0.09279831 

0.093138137 

+0.36 

Pl4 

0.00826216 

0.008079838 

-2.26 

Pl5 

0.00100513 

0.000919982 

-9.26 

P-16 

0.02671924 

0.027779444 

+3.82 

Pl7 

0.00450124 

0.004579908 

+1.72 

^18 

0.00161537 

0.001699966 

+4.98 

^19 

0.00276439 

0.002599948 

-6.32 

o 

0.19223051 

0.193556129 

+0.68 

8 

0.080323504 

0.082366904 

+2.48 

7.  CONCLUSION 

Both  simulation  studies  affirm  the  validity  of  this  analysis.  The  main  drawback  of  the  procedure  is 
that  both  its  time  and  space  complexities  are  exponential  in  M;  therefore,  it  is  impractical  for  evaluat¬ 
ing  multiplexers  with  interesting  values  of  M,  N and  B.  To  circumvent  this  difficulty,  we  propose  an 
approximation  multiplexer  model  to  evaluate  scenarios  in  which  A/,  ’s  and  M  are  large. 

8.  AN  ALTERNATE  MULTIPLEXER  MODEL 

In  the  previous  section,  we  analyzed  an  ATM  statistical  multiplexer  subject  to  heterogeneous 
groups  of  on-off  traffic  sources.  The  durations  of  the  on  and  off  periods  of  a  traffic  source  are  geo¬ 
metrically  distributed  with  parameters  a  and  /3,  respectively.  A  sample  path  of  the  state  of  a  traffic 
source  when  a  and  [3  are  close  to  1  is  illustrated  in  figure  3.  Notice  that  the  burst  periods  tend  to  last  a 
long  time;  that  is,  given  that  the  state-of-the-traffic  source  is  on  at  time  k,  then  it  is  highly  probable 
that  the  traffic  source  is  still  on  at  time  k+dk  even  if  6k  is  large.  When  cell  generation  tends  to  cluster 
together,  we  say  that  the  generation  process  is  positively  correlated  (figure  3).  The  correlativeness  of 
a  traffic  source  has  a  big  impact  on  the  cell-loss  behavior  at  a  multiplexer,  particularly  when  the 
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buffer  size  of  the  multiplexer  is  large,  since  the  large  buffer  permits  the  cells  to  interact.  Since  the  on- 
off  traffic  source  model  proposed  in  the  previous  section  can  account  for  the  correlative  nature  of  cell 
arrivals,  we  refer  to  that  on-off  model  as  the  correlated  burst  on-off  model.  Under  certain  asymptotic 
conditions,  the  correlated  burst  on-off  model  may  be  justifiably  replaced  by  a  simpler  non-correlated 
burst  on-off  model,  which  is  described  in  the  following  section. 


Figure  3.  Positively  correlated  traffic  generation. 

Consider  the  sample  paths  of  the  states  of  two  traffic  sources,  as  shown  in  figure  4.  Sample  path  #1 
is  constructed  according  to  the  correlated  burst  on-off  model  with  a  =  0.9  and  =  0.9.  Note  that 

1 


Sample  path  #2  (figure  4)  is  constructed  such  that,  at  each  time  slot,  the  traffic  source  is  on  with 
probability  0.5  and  off  with  probability  1-0.5;  that  is,  the  traffic  source  turns  on  and  off  at  some 
time,  k,  with  the  same  probability,  independent  of  the  k.  We  refer  this  on-off  traffic  model  as  the  ran¬ 
dom  burst  on-off  model.  Notice  that  over  a  long  observation  period,  the  total  durations  of  all  the  on 
periods  of  both  sample  paths  are  approximately  the  same;  however,  path  #1  (figure  4)  is  much  more 
positively  correlated,  and  therefore,  would  effect  a  much  higher  cell  loss  at  a  multiplexer. 

SAMPLE  PATH  #1 


SAMPLE  PATH  #2 


Figure  4.  Sample  paths  of  two  types  of  on-off  traffic  sources. 

Now  consider  the  superpositions  of  three  sample  paths  of  the  correlated  burst  on-off  type  and  of 
three  sample  paths  of  the  random  burst  type,  as  shown  in  figures  5  and  6.  The  two  superpositions 
appear  to  be  similar  to  each  other,  even  though  the  sample  paths  that  constituted  them  are  quite  differ¬ 
ent.  The  superposition  of  N  sample  paths  of  the  correlated  burst  on-off  model  is  stochastically  similar 
to  the  superposition  of  sample  paths  of  the  random  burst  on-off  model,  when  N  is  large.  The  correla¬ 
tiveness  of  an  individual  traffic  stream  is  smoothed  by  the  superposition.  Extensive  simulation  studies 
have  confirmed  this  interesting  phenomenon. 
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Figure  5.  Superposition  of  correlated  on-off  sources. 


Figure  6.  Superposition  of  random  on-off  sources. 


The  random  burst  on-off  traffic  streams  have  simpler  analytical  property  and,  if  they  were  used  as 
the  traffic  streams  for  a  multiplexer,  would  simplify  the  analysis  considerably;  therefore,  we  propose 
to  use  random  burst  on-off  traffic  sources  whenever  the  number  of  traffic  sources  of  a  correlated  burst 
on-off  traffic  type  is  large.  The  modified  multiplexer  model  with  random  burst  on-off  model  is  as  fol¬ 
lows: 


•  M  is  the  number  of  different  traffic  types 

•  Ni  is  the  number  of  traffic  sources  of  the  i-th  type  i=l, ...,  M 

•  Ri  is  the  number  of  cells  generated  by  an  i-th  type  traffic  source  in  a  time  slot  during  an  on 
period 

•  B  is  the  maximum  number  of  cells  that  may  be  buffered  at  the  multiplexer 

•  C  is  the  number  of  cells  that  may  be  transmitted  by  the  output  link  of  the  multiplexer  per  time 
slot;  Assume  that  1  <  C 

We  make  the  following  assumptions: 

a.  An  on  period  of  a  traffic  source  in  which  the  i-th  type  lasts  Xt  time  slots,  where  Xk  is  geometri¬ 
cally  distributed  with  parameters  at 

b.  An  off  period  of  a  traffic  source  in  which  the  i-th  type  lasts  Yt  time  slots,  where  YL  is  geometri¬ 
cally  distributed  with  parameter  1  -  a, 

c.  The  number  of  cells  at  the  buffer  at  a  time  k+1  is  computed  as  max(min(/  +  R  ■  J  -  C,B),  0), 
where  i  is  the  number  of  cells  in  the  buffer  at  time  k;  J  is  the  state  of  the  traffic  sources  at  the 
start  of  time  k 

M 

d.  /?//,«,-  <  C 
/=  1 

M 

e.  Rfii  >  C 

/=  I 

Notice  that  the  only  difference  between  the  correlated  burst  on-off  model  and  the  random  burst  model 
is  that,  in  the  correlated  model,  the  on  and  off  periods  are  determined  by  two  parameters,  a,-  and  f>i, 
whereas  in  the'random  burst  model,  the  periods  are  determined  by  only  one  parameter,  a,-. 

We  define  the  following  variables: 

•  Bk  the  number  of  cells  in  the  buffer  at  the  beginning  of  the  k-th  time  slot 

•  Jik  the  number  of  the  i-th  type  traffic  sources  that  are  on  at  the  beginning  of  k-th  time  slot 

•  S(M,£Q  =  I  0  <  j,  <  Nj  i  =  1,2,..., M) 

•  MBk+l  =  j  I  Bk  =  i,Jk  =  f)  the  probability  that  Bk+X  =  j,  given  that  Bk  =  i  and  Jk  =  j 

•  Pr{Bk+ 1  =  j  I  Bk  =  i)  the  probability  that  Bk+]  =  j,  given  that  Bk  =  i 

•  Pr(Bk+ 1  =  j)  the  probability  that  Bk+\  =  j 
We  find  that 

Pr(Bk+ 1  =  j  I  Bk  =  i,Jf  =  j)  =  |{y  =  max(min(i  +  R  ■  j_  -  C,B),0)j ,  (}1) 

Pr(Bk+i  =  j  I  Bk  =  i)=  ^  ?{y  =  max(min(i  +  R-  j_  -  C,B),  0)J  Pr(A-  =  j) ,  (12) 

jeS(M.N) 
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Pr{Bk+\  =  y)  =  j  ^  |{y  =  max(min(/  +  R  ■  j_  -  0)JPr(7t  =  ;)  j Pr(Bk  =  i)  .  (13) 

/  =  ()  \  jjES{M.l£)  J 

In  steady  state,  the  number  of  cells  in  the  buffer,  L,  is  independent  of  the  time,  k;  / ,  the  number  of  on 
traffic  sources  of  each  type  at  a  time  slot,  is  also  independent  of  the  time,  k;  therefore,  we  have  the 
following  fixed-point  relation  for  Pr(L  =  i )  for  i= 0,  1, ...,  B: 

Pr(L  =  7)  =  X(  X  £{• /  =  max(min(/  +  R  ■  j-C,B),0)\Pr(J  =  j)  )  Pr(L  =  i) ,  (14) 

where 


M  / 

PrU=j)  =  n 

The  Pr(L  =  i)’s  may  be  obtained  iteratively.  We  formulated  the  following  iteration  procedure  to  obtain 
the  distribution  function  of  the  variable,  L.  Pi  is  used  to  denote  Pr(L  -  i ). 

k  =0; 


N, 

h 


)aJ‘(l  -afrii. 


(15) 


(Pk)i  =  j  for  i=lAJ5 


( pt+')i  =  0  for 


While  ( 

{ 


PMi  ~  Pki 


5  for  0  <i  <B) 


For  0  <  j  <  B 


{ 


pick  any  zk  such  that  0  <  zk  <  1; 


pkj  =(1  -Zt)pk  lj  +EiS 


Z  qy  =  max(min(i  +  R-£  -  CJ3),0)\Pr(J  =  /) 
i€5(«yv)  [  J 


,1-1 
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s-i>v 

j= 0 

For  0  <  i  <B 


{ 


The  procedure  terminates  whenever  it  has  found  a  consistent  set  of  pkj  =  Pr(L  =j)’s. 

We  conclude  that  there  exists  a  unique  solution  to  the  iteration  scheme,  since  the  system  is  ergodic. 
Note  that  the  existence  of  a  solution  does  not  guarantee  convergence  of  the  iteration.  Extensive  com¬ 
putational  experience  suggests  that  the  procedure  always  converges.  Each  iteration  step  requires 

M 

0(B 2  ]~[  Nj)  computations;  therefore,  the  scheme  is  not  practical  for  interesting  values  of  Nt  and  M. 

i  ~  1 

We  will  propose  a  heuristic  that  would  considerably  reduce  the  computational  complexity  of  the  pro¬ 
cedure. 

Once  that  Pr(L  =  i)  for  i= 0,  1,  2, ...,  B  are  obtained,  we  may  derive  the  steady-state  cell-loss  proba¬ 
bility  of  the  multiplexer  as  follows: 

We  define  the  variables: 

•  NL(k,j)  the  number  of  cells  that  are  lost  in  a  time  slot,  given  that  there  are  k  cells  buffered  at  the 
beginning  of  the  slot  and  that  f  is  the  state  of  the  traffic  sources  at  the  beginning  of  the  time  slot 

•  NL(k)  the  number  of  cells  that  are  lost  in  a  time  slot,  given  that  there  are  k  cells  buffered  at  the 
beginning  of  the  slot 

•  NA(kj)  the  number  of  cells  that  arrive  in  a  time  slot,  given  that  there  are  k  cells  in  the  beginning 
of  the  slot  and  j_  is  the  state  of  the  traffic  sources  at  the  beginning  of  the  time  slot 

•  NA(k )  the  expected  number  of  cells  that  arrive  in  a  time  slot,  given  that  there  are  k  cells  in  the 
beginning  of  the  slot  and  j_  is  the  state  of  the  traffic  sources  at  the  beginning  of  the  time  slot 

•  NA  the  expected  number  of  cells  that  arrive  in  a  time  slot 

•  NL  the  expected  number  of  cell  losses  in  a  time  slot 

•  Cp  the  steady-state  cell-loss  probability  of  the  multiplexer 
We  find  that 


NL(kJ)  = 


K'i_~(.B  —  k+C)  B  —  k+C<R-j_ 
0  otherwise 


(16) 


NL(k)=  £  NL(k,j)  Pr(7  =  j) , 

ies(Mg) 


(17) 
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(18) 


£S 

NL  =  ^{NL(k)Pr(L  =  k)} , 

k~0 

M 

NA(k,j)  =  R  •  i  =  YR'J'> 

i~  1 

NA(k)  =  Y  NA(k,j)Pr(J  =/)  =  2  (<&  •  j)PrU  =  ;)} , 

NA  =  X  {(£-Z>Pr^  =  j)}- 


The  cell-loss  probability  during  a  time  slot  is 


Cl  ~ 


B  1 

I 

f  V  1 

X  (R~  i~  (B  -  k  +  C))Pr(J=M(R-j  ZB-k+C)\ 

\Pr(L  =  k) 

NL  _  *-°  1 

[ye5(M.AD  1 

NA 


(19) 

(20) 
(21) 


(22) 


X  {(R  ■  j)Pr(J  =  j)} 

;g5(m.ao 

From  the  above  expression,  the  cell-loss  probability  may  be  obtained  in  principle;  however,  the  com¬ 
putational  complexity  is  exponential  to  the  sizes  of  the  parameters  A'-’s,  Ri’s,  and  M.  We  propose  a 
heuristic  that  will  make  the  requisite  computation  feasible. 

Notice  that 


B  I  \  B 

X  X  =  +  R  'l  -  C,5),0)]Pr(i  =  j)  )  =  X  X  PrU  =  &> 


(23) 


i  =  0  \  j£S(M,N) 


where 


Sijj  -  {i  I  j  =  max(min(7  +  R-  [-  C,B),0)} . 

For  a  fixed  pair  the  set,  SUi,  should  be  much  smaller  than  S(M,N).  We  develop  a  procedure  that 
enumerates  the  set  SUi  without  testing  each  j  E  S(M,N). 

We  partition  SUi  into  SUj\J  =  0],  Sufil  <  j  <  B]  and  S.JB  <  j],  where 


SijjU  =  0 ]  =  {;£  S(M,N)  :  R  ■  j_  <  C-i) , 

(24) 

Syjl  =  {I  e  WAD  :  R  ■  j  =  C  +  j-i] , 

(25) 

S.JB  =  j]  =  (/G  WAD  '■  R  '  l~  C  +  B-i } . 

(26) 

We  define  the  following  sets: 


A  ] 

*  ^fl...Af][—  &]  “  ' 

j  G  S(A/,AQ 

• 

1! 

II 

i  6  S(A/,tf) 

M  1 

X^'/  =  * 

j)  £  S(M,N)  1 

Af  1 
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We  formulate  a  procedure  that  will  enumerate  the  sets  k]  and  s[UJfl[;>  i]  without 

testing  each  j  E  S(M,N). 

The  notations  are  as  follows: 


•  S/[=  k]  =  {y  I  0  <y,  <  N,,Rj  =  k], 

•  A  X  B  =  {( ij )  I  i  E  AJ  E  B) , 


/^,  X  0'})  M  >  m  >  1 

“  jesf=k) 

{y  I  0  <  y  <  Nh  jR[  =  k}  m=  1 


The  summation  symbol  represents  the  wn/orc  set  operation. 


We  define 


•  S,[<  k]  =  (y  I  0  <y  <  N,t  Rj  <  k} 


•  *] 


X  ( k-jRm )]  X  {y})  m  >  m  >  1 

‘  JSSjlsk) 

{j\0<j<NhjRt<k}  m=  1 


We  introduce  the  notation 


•  Sj  =  {/  I  0  <  y  <  A?,} , 


^(•5[i..jn-il[—  (M?«M  X  {y}) 
jesm 


[{/  I  0  <  j  <  N>,  jRt  >  k} 


M  >  m  >  1 
m  =  1 


Using  this  enumeration  scheme,  the  summation  Of  the  iterative  procedure  may  be  executed 

l£S(,M.N) 

without  testing  each  y  E  S(M,N).  This  reduces  computation  considerably.  In  the  worst  case,  our  com- 

M 

putational  procedure  to  calculate  the  exact  cell-loss  probability  is  still  0(B2]~]yv,),  although  its  aver- 

/'=  1 

age  complexity  is  much  lower.  For  the  case  that  Nfs  are  large,  we  develop  the  following  efficient  pro¬ 
cedure  to  approximate  the  cell-loss  probability  of  the  multiplexer. 


9.  DIFFUSION  APPROXIMATION 


When  Ni  is  large  for  each  i=l,  2, ...,  M,  diffusion  approximation  may  be  used  to  estimate  Pr(J  =  j). 
If  Ni  is  large,  Np. ,  §>  1  and  Np.;(  1  —  a,)  S>  1,  then  by  the  Laplace-DeMoivre  Theorem, 


Pr(J,  =  y)  ~ 


fhio) 


-U-M/r 


2a? 


where  p,  =  and  of  =  N,a,{\  -  a ,). 


(27) 


17 


We  approximate  the  discrete  random  variable,  Jt,  with  a  continuous  variable,  The  probability 
distribution  for  Yi  is 


Pr(Y,  <  x)  = 


where  K ,  = 


0  <  x  <  N  i , 


(28) 


Yi  is  approximately  Gaussian. 
Finally,  if  we  define 


(29) 

/  =  ] 

M  M 

Y is  also  approximately  Gaussian  with  a  mean  of  fi  -  X  Rfip.,  and  a  variance  o2  =  X  -  a,). 

i= 1  i=l 

The  distribution  function  of  Y  is 


Pr(Y 


X 

~x)  =  wfepj' 


-(T-fl)2 


2o2 


M 

dr  0<*<X  , 


/=  1 

where  £  =  — -1  , 

Jlno2  I  ‘ 


2  a2 


dr . 


(30) 


We  conclude  that 


X  Pr(2  =  J)  ~  =  OJ/VfF  ^  C  -  0 

iesuj 

+  £{0  <j  <  B }Pr(C  +  j-i-l  <  Y  <  C  +  j-i)  +  £{B  =  j}(l-Pr(Y  <  C  +  B-i)) . 


(31) 


The  Pr(L  =  i)’s  may  be  obtained  from  the  following  iteration  scheme  in  which  ^  P(J  =  j )  is 

iesiJi 

approximated  by  the  expression  of  equation  (31). 
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k  =  0; 


(Pk)i  =  j  for  i'=U,..3 
=  0  for  1=1,2,. .,3 
While  ( 


p*+\-  -  pki  \  >  8  for  0  <  i  <  B) 


k  =  k+l  ; 

For  0  £  j  <B 

{ 

pick  any  £k  such  that  0  <tt  <  1 ; 

Pkj  =  (1  "  lk)pk~'j  +  e*  2  Pr{J_=i)pk-\ 


*«2>V 

M 

For  0  <•  i  <  B 


{ 


We  also  may  approximate  NL  and  NA  as  follows: 

We  note  that 

X  {( E.'j-(B-k  +  C))Pr(J  =j)}=  y\{(R-j-(B-k+  C))Pr(J  =  j 

jesw.®  jBS  . 

-  kj 

where  Sk{  m  {/  e  S(M,N)  :B-  k  +  C  <  R-j\. 


(32) 
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Again,  we  may  approximate  the  sum 


1V1 


X  {(£  '  HB-k  +  O)  Pr (J  =  ;)}  ~  [ 

ieskj  J 


( *-ak)fJj)dT , 


where  ak  -  B  -  k  +  C  and  f/r) 
Finally, 

The  cell-loss  probability  is 


Kjljto2 


-{r-fxY 


2*2 


m 

NA  ~  . 

f=l 

M 

JjRiNi 

B  1 

£(  (  (r-ak)fyir)dr)Pv(L  =  k) 

*=()  I 


CL~ 


(33) 


(34) 


(35) 


/=1 

Using  this  approximation,  the  computation  at  each  iteration  step  is  reduced  to  the  evaluation  of  B 2 

M 

error  function  integrals.  The  complexity  of  evaluating  an  integral  is  G>(^T  /?//,.);  therefore,  the  overall 

f=l 

M 

complexity  of  the  approximation  procedure  is  0(B 2  ^  R,N,).  If  one  evaluates  each  integral  using  a 

/=  1 
M 

fixed  number,  rather  than  C>(]T  /e,A^,),  of  steps,  the  complexity  of  the  algorithm  is  reduced  to  0(B2).  We 

/'=  1 

performed  several  experiments  using  the  0(B2)  algorithm  and  found  that  it  yielded  virtually  the  same 

M 

results  as  the  0(B 2  Rfl,)  algorithm. 


10.  NUMERICAL  EXAMPLES  II 

We  conducted  several  simulation  studies  to  evaluate  the  accuracy  of  exact  and  approximation  algo¬ 
rithms  that  estimate  the  cell-loss  probability  of  the  random-burst  multiplexer  model.  The  studies 
showed  good  agreement.  Tables  3  and  4  list  the  results  from  two  of  the  simulation  studies. 

10.1  EXAMPLE  #1  (EXACT  ANALYSIS) 

B  =  10 
C  =  50 
M  =  2 
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Nj  =  3  aj  =  0.2  Rj  =  40 
N2  =  10  a.2  =  0.8  R2  =  1 


We  denote  P,  as  the  steady-state  probability  that  the  buffer  content  of  the  multiplexer  at  the  start 

of  a  time  slot  is  i  and  e  as  the  steady-state  cell-loss  probability.  Relative  error  is  defined  as 

,  _  simulation  -  analysis 

100  x - ; — — - - — . 

simulation 


Table  3.  Exact  algorithm. 


Analysis 

Simulation 

Relative  Error  (%) 

Po 

0.82978 

0.83290 

-0.376 

Pt 

0.00106 

0.00110 

-3.78 

P2 

0.00148 

0.00150 

-1.35 

Pz 

0.00228 

0.00220 

+3.51 

Pa 

0.00350 

0.00330 

+5.71 

Ps 

0.00524 

0.00500 

+4.58 

Pe 

0.00816 

0.00780 

+4.41 

P7 

0.01230 

0.01170 

+4.88 

Ps 

0.01444 

0.01440 

+0.28 

Ps 

0.01150 

0.01770 

-1.74 

P10 

0.110260 

0.10850 

+1.60 

8 

0.1069174 

0.105817096 

+1.03 

10.2  EXAMPLE  #2  (APPROXIMATION  ANALYSIS) 

B  =  20 
C  =  3556 
M  =  5 

Nj  =  50  Rj  =  10  aj  =0.6 
N2  =  50  R2  =  30  a2  =  0.7 
N3  =50R3  =  20a3  =  0.5 
N4  =  50  R4  =  50  a.4  =  0.3 
N5  =  50R5  =  60a5  =  0.2 


Table  4.  Approximation  algorithm. 


Analysis 

Simulation 

Relative  Error  (%) 

Po 

0.9096 

0.9079 

+0.19 

Pa 

0.0050 

0.0060 

-20.00 

Pu 

0.0060 

0.0057 

+5.00 

P20 

0.0794 

0.0802 

-1.01 

e 

0.003015543 

0.002997045 

+0.61 
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11.  APPLICATION:  ATM  EFFICIENCY 

How  efficient  is  ATM?  This  is  one  of  the  most  frequently  asked  questions  by  ATM  network 
designers.  The  following  example  illustrates  how  the  algorithms  formulated  in  the  preceding  sections 
would  aid  in  addressing  this  question. 

To  simplify  our  discussion,  we  consider  an  ATM  network  consisting  of  only  one  transmission  link. 
The  network  is  to  support  many  connections,  each  carrying  the  same  traffic  type.  We  wish  to  deter¬ 
mine  what  is  the  maximum  link  efficiency  achievable  such  that  the  quality  of  service  requirement 
(QOS)  of  each  connection  supported  is  satisfied.  The  answer  depends  on  the  following  parameters: 

•  C  the  capacity  of  the  transmission  link 

•  B  the  size  of  the  buffer 

•  N  the  number  of  connections  being  supported 

•  TP  the  traffic  profile  of  the  connections 

•  Q  the  QOS  requirement  of  the  connections 
We  make  the  following  simplifying  assumptions: 

•  The  traffic  profile,  Z£  is  specified  in  terms  of  the  peak  traffic  rate,  the  mean  traffic  rate,  and 
the  expected  burst  period  of  a  connection. 

•  The  QOS  requirement,  Q,  is  usually  defined  in  terms  of  both  the  maximum  cell-loss  rate  and 
the  maximum  end-to-end  delay  tolerable  by  the  connection.  To  simplify  the  discussion,  we 
assume  that  the  QOS  requirement  is  measured  only  in  terms  of  the  maximum  cell-loss  rate 
tolerable  by  a  connection  . 

The  transmission  link  is  statistically  shared  ( multiplexed )  among  N  connections.  When  the  total 
cell-generation  rate  of  the  N  connections  is  below  the  capacity  of  the  link,  C,  no  loss  occurs;  however, 
when  the  instantaneous  total  traffic  rate  exceeds  C,  the  excess  cells  will  be  stored  in  the  buffer.  When 
the  buffer  is  full,  cells  will  be  dropped. 

We  will  use  the  following  examples  to  illustrate  how  the  parameters  influence  the  answer  to  our 
question. 

For  each  of  the  two  examples  below,  we  assume  the  following: 

The  network  is  to  support  packetized  voice  connections.  We  assume  that  the  voice  encoder  of  a  con¬ 
nection  would  sample  the  incoming  voice  signal  at  8000  samples/sec  and  that  each  sample  is  encoded 
into  an  8-bit  unit.  A  cell  is  packed  with  48  samples  (bytes).  We  also  assume  that  the  encoder  has  a 
silence  detector  and,  therefore,  would  not  generate  samples  during  the  period  that  the  talker  is  silent. 
We  assume  that  the  expected  burst  period,  the  interval  during  which  the  talker  speaks,  is  0.96  sec-  . 
onds;  the  expected  silence  period,  the  interval  during  which  the  talker  is  silent,  is  1.60  seconds. 

Based  on  the  above  assumptions,  the  traffic  profile  of  each  voice  connection  is  as  follows: 

•  The  peak  cell  rate  of  the  connection  is  167  cells/sec. 

•  The  mean  cell  rate  of  the  connection  is  62  cells/sec. 

•  The  expected  burst  period,  the  period  during  which  the  voice  source  would  generate  traffic  at 
its  peak  rate,  is  0.96  seconds. 
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•  It  is  assumed  that  the  cell-loss  requirement  of  each  connection  is  10‘7  and  that  an  overall  cell- 
loss  rate  of  10'8  at  the  transmission  link  suffices  to  effect  a  cell-loss  rate  of  10‘7  to  each  connec¬ 
tion. 


11.1  EXAMPLE  1 

We  assume  that  the  link  can  transmit  up  to  10  Mbits/sec  (-26000  cells  per  second). 

The  buffer  may  store  up  to  100  cells. 

What  is  the  maximum  efficiency  achievable  by  the  link?  That  is,  what  is  the  maximum  number  of 
voice  connections  that  can  be  supported  by  the  link  such  that  the  cell-loss  requirement  of  each  con¬ 
nection  is  satisfied? 

We  use  the  multiplexer  model  with  correlated  burst  on-off  traffic  sources.  According  to  the  model, 
the  maximum  number  of  connections  that  can  be  supported  by  the  link  is  N=  310.  The  efficiency  is 


11.2  EXAMPLE  2 

We  retain  the  assumptions  made  in  example  1,  with  the  exception  that  the  link  capacity  is  only 
1  Mbits/sec. 

The  answer  to  the  same  question  is  N  =  20.  The  efficiency  is  ~  —  x  —  =  0.48. 

2600 

Notice  that  the  network  in  example  1  has  a  much  higher  statistical  gain  than  that  of  example  2. 
This  reenforces  the  proven  principle  that  networks  with  higher  capacity  transmission  links  can 
achieve  significantly  higher  statistical  gain  than  low-capacity  networks. 

12.  FUTURE  WORK  AND  SUMMARY 

In  this  report,  a  discrete-time  ATM  statistical  multiplexer  with  correlated  burst  on-off  traffic 
sources  is  proposed.  Algorithms  to  compute  both  the  exact  steady  state  and  transient  cell-loss  proba¬ 
bilities  of  the  multiplexer  have  been  formulated.  Simulation  studies  confirm  the  validity  of  the  exact 
algorithms.  Unfortunately,  the  computational  requirement  of  the  algorithms  prohibits  their  use  in 
evaluating  all  ATM  multiplexing  scenarios.  As  a  remedy,  we  proposed  a  multiplexer  model  with  ran¬ 
dom  burst  on-off  traffic  sources.  This  model  has  similar  cell-loss  behaviors  as  the  correlated  burst 
model’s  when  the  number  of  traffic  sources  of  each  traffic  type  being  multiplexed  is  large.  For  the 
random  burst  model,  we  formulated  an  algorithm  to  compute  the  multiplexer’s  cell-loss  rate  (transient 
and  steady  state).  The  time  complexity  of  the  exact  algorithm  is  still  exponential;  for  this  reason,  we 
developed  an  efficient  approximation  algorithm  to  compute  the  cell-loss  rate.  Extensive  simulation 
studies  had  been  performed  to  evaluate  the  accuracy  of  the  approximation  algorithm,  and  for  all 
parameters  of  practical  interest,  the  algorithmic  results  and  the  simulation  outputs  showed  perfect 
agreement.  Finally,  we  illustrated  the  practical  use  of  the  algorithms  presented  in  this  report  to  deter¬ 
mine  the  expected  efficiency  of  an  idealized  ATM  network. 

The  multiplexer  model  with  correlated  burst  on-off  traffic  sources  is  considered  to  be  one  of  the 
most  general  multiplexer  models  that  has  ever  been  proposed.  As  far  as  we  are  able  to  determine,  no 
exact  solutions  (closed-form  or  algorithmic)  to  comparably  general  models  are  known.  The 
contribution  of  this  paper  is  the  development  of  exact,  easily  implementable,  and  numerically  stable 
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algorithms  to  compute  the  transient  and  steady-state  cell-loss  probabilities  of  a  faithful  ATM  statisti¬ 
cal  multiplexer  model.  The  main  drawback  of  the  solution  is  that  it  is  algorithmic  rather  than  in 
closed-form.  An  obvious  future  work  on  this  subject  is  to  find  a  closed-form  solution  for  the  cell-loss 
probability  of  this  model. 
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